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ABSTRACT 

. The sources and fluxes of superGZK neutrinos, E > 10 20 eV, are discussed. The 

| fluxes of cosmogenic neutrinos, i.e. those produced by ultra-high energy cosmic 

£N| . rays (UHECR) interacting with CMB photons, are calculated in the models, 

which give the good fit to the observed flux of UHECR. The best fit given in 
flY 1 . no-evolutionary model with maximum acceleration energy E ma ^ = 1 x 10 21 eV 

^0 ' results in very low flux of superGZK neutrinos an order of magnitude lower 

. than the observed flux of UHECR. The predicted neutrino flux becomes larger 

fN| ' and observable by next generation detectors at energies 10 20 — 10 21 eV in the 

evolutionary models with E max = 1 x 10 23 eV. The largest cosmogenic neutrino 
flux is given in models with very flat generation spectrum, e.g. oc E~ 2 . The 
neutrino energies are naturally high in the models of superheavy dark matter 
' and topological defects. Their fluxes can also be higher than those of cosmogenic 

\Q , neutrinos. The largest fluxes are given by mirror neutrinos, oscillating into 

0^ ' ordinary neutrinos. Their fluxes obey some theoretical upper limit which is very 

weak, and in practice these fluxes are most efficiently limited now by observations 
' of radio emission from neutrino-induced showers. 
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1. Introduction 



The abbreviation 'SuperGZK neutrinos' implies neutrinos with energies above 
the Greisen-Zatsepin-Kuzmin d cutoff E GZ k ~ 5 x 10 19 eV. Soon after theoretical 
discovery of the GZK cutoff, it has been realized that this phenomenon is accompanied 
by a flux of UHE neutrinos, which in some models can be very large In 80s it was 
understood that topological defects can produce unstable superheavy particles with 
masses up to the GUT scale ^ and neutrinos with tremendous energies can emerge 
due to this process^. 

It has been proposed that SuperGZK neutrinos can be detected observing the 
horizontal Extensive Air Showers (EAS) The exciting prospects for detection 
of SuperGZK neutrinos have appeared with the ideas of space detection, e.g. in 
the projects EUSOE* and OWL . The basic idea of detection can be explained by 
example of EUSO. 

The superGZK neutrino entering the Earth atmosphere in near-horizontal direc- 
tion produces an EAS. The known fraction of its energy, which reaches 90% , is 
radiated in form of isotropic fluorescent light. An optical telescope from a space 
observatory detects this light. Since the observatory is located at very large height 
(~ 400 km) in comparison with thickness of the atmosphere, the fraction of detected 
flux is known, and thus this is the calorimetric experiment (absorption of light in 



the upward direction is small). A telescope with diameter 2.5 m controls the area 
~ 10 5 km 2 and has a threshold for EAS detection E th ~ 1 x 10 20 eV. 

The very efficient method of superGZK neutrino detection is given by observations 
of radio emission by neutrino-induced showers in ice, salt and lunar regolith. This 
method has been originally suggested by G. Askaryan in 60s Propagating in the 
matter the shower acquires excessive negative electric charge due to involvement of 
the matter electrons in knock-on process. The coherent Cerenkov radiation of these 
electrons produces the radio pulse. Recently this method has been confirmed in 
the laboratory measurements ^ . There were several searches for such radiation from 
neutrino-induced showers in the Antarctic and Greenland ice and in the lunar regolith. 
In all cases the radio-emission can be observed only for neutrinos of extremely high 
energies. The upper limits on the flux of these neutrinos have been obtained: in 
GLUE experiment by radiation from the moon, in FORTE experiment by 
radiation from the Greenland ice and in RICE experiment from the Antarctic ice. 

The characteristic feature of the detection methods described above is the high 
energy threshold, typically E > 1 x 10 19 — 1 x 10 20 eV. How neutrinos of these energies 
can be produced? 

The most conservative mechanism of superGZK neutrino production is py mech- 
anism of collisions of accelerated protons/nuclei with low-energy CMB photons. To 
provide neutrinos with energies higher that 1 x 10 20 eV the accelerated protons must 
have energies higher (or much higher) than 2 x 10 21 eV. For shock acceleration this 
energy can reach optimistically 1 x 10 21 eV. One has raise his hopes on less devel- 
oped ideas of acceleration such as acceleration in strong e-m waves, exotic plasma 
mechanisms of acceleration and unipolar induction. 

The top-down scenarios can easily provide neutrinos with energies higher and 
much higher than 1 x 10 20 eV. The idea common for many mechanisms is given by 
existence of superheavy particles with very large masses up to GUT scale. In Grand 
Unified Theories (GUT) these particles (gauge bosons and higgses) are short-lived. 
In the cosmic space they are produced by Topological Defects (TDs). The decay of 
these particles results in the parton cascade, which is terminated by production of 
pions and other hadrons. Neutrinos are produced in their decays. 

The superheavy particles are naturally produced at post-inflationary stage of the 
universe. The most reliable mechanism of production is gravitational one. The masses 
of such particles can reach 10 13 — 10 14 GeV. Protecting by some symmetry (e.g. 
gauge symmetry or discrete gauge symmetry like R-parity in supersymmetry), these 
particles can survive until present cosmological epoch and produce neutrinos in the 
decays or annihilation. 



2. Upper limits on superGZK neutrino flux 



There are two different upper limits on UHE neutrino fluxes: cascade upper 



limit ^ and cosmic ray upper limits (Waxman-Bahcall ■ ■ and Mannheim-Protheroe- 
Rachen^J). The cosmic ray upper limits are not relevant for superGZK neutrinos 
because this limit is not valid for top-down scenarios and it is automatically satisfied 
for cosmogenic neutrinos, since their fluxes are calculated in the models which explain 
the observed UHECR. 

The cascade upper limit on HE and UHE neutrino fluxes EE3 j s provided due to 
e-m cascades initiated by HE photons or electrons which always accompany produc- 
tion of neutrinos. Colliding with the target photons, a primary photon or electron 
produce e-m cascade due to reactions 7 + 7 tar — > e + + e~, e + 7 tar — > e' + 7', etc (see 
Fig. Q). The standard case is given by production of HE neutrinos in extragalactic 
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Figure 1: Developing of e-m cascade in collisions with background target (t) photons. 

space, where cascade develops due to collisions with CMB photons (7 tar = 7cmb)- In 
case the neutrino production occurs in a galaxy, the accompanying photon can either 
freely escapes from a galaxy and produce cascade in extragalactic space, or produce 
cascade on the background radiation (e.g. infra-red) within the galaxy. In the latter 
case the galaxy should be transparent for the cascade photons in the range 10 MeV 
- 100 GeV. 

The spectrum of the cascade photons is calculated I 5 | 15 | l0 } : j n i QW ener gy p ar t 
it is oc E~ 3 / 2 , at high energies oc E~ 2 with a cutoff at some energy e 7 . The energy 
of transition between two regimes is given approximately by e c ~ (e t /3)(e 7 /m e ) 2 , 
where e t is the mean energy of the target photon. In case the cascade develops 
in extragalactic space e t = 6.35 x 10~ 4 eV, e 7 ~ 100 GeV (absorption on optical 
radiation), and e c ~ 8 MeV. The cascade spectrum is very close to the EGRET 
observations in the range 3 MeV - 100 GeV 13. The observed energy density in this 
range is c^egret ~ (2 — 3) x 10~ 6 eV/cm 3 . The upper limit on HE neutrino flux 
J u (> E) is given by chain of the following inequalities 

An r°° An r°° Atx 

u cas > — / EJ u (E)dE > —E / J u {E)dE = —EJ U {> E), 
c Je c Je c 

which in terms of the differential neutrino spectrum J U (E) results in 



E 2 J V (E) < — u cas , with ^ cas < Regret- 



(1) 



Unless otherwise is stated, here and everywhere below the neutrino flux J u is given 
as sum of all neutrino flavors. 

Eq. gives the rigorous upper limit on the neutrino flux. It is valid for neutrino 
production by HE protons, by TDs, by annihilation and decays of superheavy parti- 
cles, i.e. in all cases when neutrinos are produced through decay of pions and kaons. 
It is valid for production of neutrinos in extragalactic space and in galaxies, if they 
are transparent for the cascade photons. It holds for arbitrary neutrino spectrum 
falling down with energy. If one assumes some specific shape of neutrino spectrum, 
the cascade limit becomes stronger. For example, for E~ 2 neutrino spectrum one 
immediately obtains 

E 2 J V {E) < f . U ~ (2) 
47r \n(E max / E min ) 

3. Cosmogenic neutrinos 

Cosmogenic has a meaning "produced by cosmic rays" . 

The most efficient mechanism for production of cosmogenic superGZK neutrinos is 
given by py collisions of protons with CMB photons: for E u > 1 x 10 20 eV the energy 
of the parent protons E p ~ 20E U is enough for photopion production in collisions with 
CMB photons. The space density of CMB photons (412 cm -3 ) is usually much larger 
than number density of the gas and optical/IR photons in the sources and outside. 

We shall reproduce here the historically first calculations^ of the diffuse neutrino 
flux produced by UHE protons colliding with CMB photons. 

Consider the universe filled uniformly by UHECR sources with space density n s 
and UHE proton luminosity L p . We assume the cosmo logical evolution of the sources 
in the form C = Cq(1 + z) m , where Co = L p n s is emissivity at the epoch with redshift 
z = and factor (1 + z) m describes the evolution. 

The production rate of a source is given as 

Q g en = (7 9 -2)L P £-^, (3) 

where all energies are measured in GeV and E min is assumed to be ~ 1 GeV. 

The diffuse UHE proton flux can be calculated from particle conservation, assum- 
ing the generation energy E g = E g (E, z) due to energy losses and integrating over all 
epochs of generation: 

J P (E) = ^-( 7s -2)£ / dt(l + z) m Ep(E,z)dE g /dE. (4) 

Unmodified diffuse spectrum is calculated with only adiabatic energy losses in- 
cluded: E g (E, z) — (1 + z)E and dE g /dE — 1 + z. Using the connection of cosmo- 
logical time t and redshift z as dz/dt = H (l + z)X(z), where 



x(z) = ^J(i + z) 3 n m + (i + z) 2 n r + n A , 



(5) 



Hq is the Hubble constant and Q m , Q r and Q\ are cosmological density in units 
of critical density of non-relativistic dark matter m, relativistic dark matter r and 
that due to vacuum energy density A, respectively one obtains for the unmodified 
spectrum of UHE protons: 

c C 

Jumn(E) = — (7 g - 2)-^-E" y9 ri ev (m, ^ max ), (6) 
where r/ ev (m, z max ,7g) is the evolutionary factor given by 

V«,(™,z™*,lg)= / dz(l + z) m -^/X(z). (7) 

It is easy to express the neutrino diffuse flux through unmodified proton flux, assuming 
that a proton undergoes several collisions with CMB photons: 



where 2/3 accounts for probability of charged pion production, 3 - for 3 neutrinos 
produced in the chain of pion decay, E u /E p fa 0.05 is a fraction of proton energy 
transferred to neutrino, and a is a fraction of energy lost by the proton in jry collision 
(a varies from 0.22 in A-resonance to 0.5 at extremely high energies); the term with 
a in Eq. (jHJ) describes approximately the subsequent jry collisions of UHE proton. 

The low-energy edge of neutrino spectrum (jSJ) is determined by energy of a proton 
at epoch z, at which the energy loss due to pion production becomes less than that 
due e + e~-pair production. 

Neutrino flux given by Eq. (jHJ) strongly depends on the parameters of cosmological 
evolution of the sources, m and £ max , as it is seen from Eqs. © and (J7J). 

The accuracy of neutrino-flux calculations by method of Ref. ® can be compared 
with exact calculations of Ref. ) where all details of py interaction were included 
and fluxes were computed for all neutrino flavors separately. The total neutrino 
moments (yields) Z pl {^ g ) calculated in Ref. must coincide with the coefficient in 
front of J unm (E) in the rhs of Eq. (JHJ). For 7^ < 2.7 the agreement is indeed 20 - 40 
%. 

The detailed calculations of UHE neutrino fluxes in the similar models with evo- 
lution of the sources and with normalization of the fluxes by the proton component 
have been performed in Refs. - ED. In this paper I will present fluxes calculated 

in the BGG models E21 ; which describe precisely the observed UHECR spectra, 
with the dip at 1 x 10 18 — 4 x 10 19 eV as the most prominent feature. 

In fact there are several BGG models which give good fit to UHECR spectra 
as measured by AGASA, HiRes, Fly's Eye and Yakutsk detectors. They differ by 
cosmological evolution of the sources, described by factor (1 + z) m , by exponent 7 9 
of the generation spectrum and by maximum energy of accelerated protons E max . 



The model with the best fit (see 1 ^) corresponds to m = (the non-evolutionary 
model), with the generation index 7 g = 2.7 and -E max = 1 x 10 21 eV. Less important 
assumption is flattening of generation spectrum at E < E c , with E c ~ 1 x 10 18 eV, 
which is needed to describe correctly the mass composition observed at E < 1 X 
10 17 eV^*. The spectrum of CR in this model is shown in Fig. |21 (upper panel) in 
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Figure 2: UHE neu trin o fluxes ^3 in the non-evolutionary (upper panel) and evolutionary (lower 
panel) BGG models E3. In the upper panel the neutrino flux accompanying the observed UHECR 
flux is shown by curve Ei^. The following parameters are used in calculations: m = 0, z ma x — 
2, 7s = 2.7, E max = 1 x I0 21 eV, E c = 1 x I0 18 eV and emissivity C = 3.5 x 10 46 erg/Mpc 3 yr. 
In the low panel the neutrino flux is maximized by the following choice of parameters: m = 4.0, 
z max = 6.0, j g = 2.45 and emissivity Cq = 1.2 x 10 46 erg Mpc^ 1 yr~ 1 . Note that fit to observational 
data in case of the evolutionary model is worse than for non-evolutionary model. 



comparison with Akeno- AGASA data. The emissivity of the sources needed to fit the 
the observed flux is Co = 3.5 x 10 46 erg/Mpc 3 yr, which corresponds to luminosity 
of a source L p = 3.7 x 10 43 erg/s for space density of the sources (powerful AGN) 



n s = 3 x 10 5 Mpc 3 . One can notice the precise agreement with the Akeno-AGASA 
data at 1 x 10 18 -8 x 10 19 eV. The explanation of the AGASA excess at E > 1 x 10 20 eV 
needs the additional CR component of another origin. The calculated neutrino flux 
is shown by curve Sz/j for sum of all neutrino flavor. This is the lowest neutrino flux 
compatible with the observed UHECR flux, because including evolution and increasing 
-Emax one increases the neutrino flux. The predicted flux of superGZK neutrinos at 
E > 1 x 10 20 eV is hardly detectable by the methods discussed above. 

Thus, the observed UHECR flux does not guarantee the detectable flux of su- 
perGZK neutrinos. 

In Fig. 121 (lower panel) neutrino flux is maximized for the BGG models, using the 
cosmological evolution, compatible with the observed UHECR flux, namely C(z) = 
(1 + z) m C at z < z max with m = 4, z max = 6 and £ = 1.2 x 10 46 erg/Mpc 3 yr. 
Very large E max = 1 x 10 23 eV is used to increase further neutrino flux. Note, that 
though the local emissivity at z = is relatively low, it was much larger in the past. 
The maximum acceleration energy E ma , x cannot be also considered as realistic, and 
used mostly for illustration of range in the predicted fluxes. Finally, this model fits 
observed UHECR spectrum worse than non-evolutionary model and has the problems 
with observed mass composition of cosmic rays at E < 1 x 10 18 eV. 

In the evolutionary models with large m and z max , which explain the observed 
UHECR, the fluxes of cosmogenic neutrinos are observable by future detectors. 

27,5 

27,0 

26,5 
"C 26,0 
" E 25,5 



==> 24,5 
"uj 

g>°24,0 
23,5 
23,0 

16 17 18 19 20 21 22 23 

log io (E/eV) 

Figure 3: UHE neutrino fluxes in the evolutionary and non-evolutionary models with proton 
generation spectrum oc E~ 2 and maximum acceleration energy E mayi = 1 x 10 23 eV. The calculated 
proton spectra (solid curves with dots) for m = and m = 4 are shown in comparison with the 
Akeno-AGASA data. The figures on spectrum curves (2, 4 and 6) show z max . The emissivities are 
Cq = 2.65 x 10 45 erg Mpc _1 yr _1 for m — and Cq = 2.2 x 10 45 erg Mpc _1 yr _1 for m = 4. 




The largest cosmogenic neutrino flux can be obtained in the models with flat 
proton generation spectra. This class of models cannot describe the spectrum at 
energies 1 x 10 18 — 1 x 10 19 eV and thus corresponds to transition from galactic to 
extragalactic cosmic rays at the ankle E ~ 1 x 10 19 eV. It has been proposed 
that large neutrino flux can be a signature of the ankle as transition from galactic to 
extragalactic cosmic rays. The calculated neutrino spectra for the proton generation 
spectrum oc E~ 2 and -Emax = 1 x 10 23 eV is presented in Fig. OJfor non-evolutionary 
model m = and for the evolutionary models with with m = 4 and z max = 2, 4 and 
6. The emissivity Cq at z = varies from 2.2 x 10 45 erg Mpc~ 3 yr _1 for m = 4 to 
2.7 x 10 45 erg Mpc~ 3 yr _1 for m = 0. 

The largest neutrino flux in Fig. E] (m — 4 and 2 max = 6) almost saturates the 
cascade upper limit: c<j cas = 1.5 x 1CT 6 eV/cm 3 , to be compared with the EGRET 
upper limit u; cas ~ 2 x 10~ 6 eV/cm 3 . 

Fig. El illustrates the range of predictions for UHE neutrino fluxes in case of flat 
generation spectrum: from modest flux in the non-evolutionary case m = to flux 
30 times higher in case of evolution. However, in all cases the maximum acceleration 
energy E m£LX = 1 x 10 23 eV is at least two orders above that obtained in realistic 
models. The top-down scenarios considered in the next sections provide these high 
energies naturally. 

4. UHE neutrinos from Superheavy Dark Matter (SHDM) 

SHDM is one of the models for cosmological cold dark matter l ^ | 28 ) rp^g mos ^ 
attractive mechanism of production is given by creati on of superheavy particles in 
time- varying gravitational field in post-inflation epoch 1 ^9 1 30 1 Creation occurs when 
the Hubble parameter is of order of particle mass H(t) ~ mx- Since the maximum 
value of the Hubble parameter is limited by the mass of the inflaton H(t) < m^ ~ 
10 13 GeV, the mass of X-particle is limited by m^, too. For example, mx ~ 3 x 
10 13 GeV results in Qxh 2 ~ 0.1, as required by WMAP measurements. 

Being protected by some symmetry, SHDM particles with such masses can be 
stable or quasi-stable. In case of gauge symmetry they are stable, in case of gauge 
discrete symmetry they can be stable or quasi-stable. Decay can be provided by 
superweak effects: wormholes, instantons, high- dimension operators etc. 

Like any other form of cold dark matter, X-particles are accumulated in the halo 
with overdensity 2.1 x 10 5 . 

SHDM particles can produce UHECR and high energy neutrinos at the decay 
of X-particles (when the protecting symmetry is broken) and at their annihilation, 
w hen the sy mmetry is exact. The scenario with decaying X-particles was first studied 
in l ^7 | 31 | 32 | ^ n interesting scenario with stable X-particles, when UHE particles 
are produced by annihilation of X-particles has been put forward in In this 

scenario superheavy X-particles have the gauge charge and they are produced at post- 



inflationary epoch by close pairs, forming the bound systems. Loosing the angular 
momentum, these particles inevitably annihilate in a close pair. 

The UHE particles (protons, pions and neutrinos from the chain of pion decays) 
are produced as a result of QCD cascading of partons. The calculations of fluxes and 

FTTl 

spectra are nowadays reliably performed by Monte Carlo and using the DGLAP 



equations 



E3 _EHJ 



The spectra of protons, photons and neutrinos are shown in Fig. H] 



for the case of SHDM particles with mass Mx = 1 x 10 14 GeV. One can observe the 
large fluxes of superGZK neutrinos with very high energies in excess of 10 22 eV. 
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Figure 4: Spectra of neutrinos (upper curve), photons (middle curve) and protons (two lower curves) 
in SHDM model compared with AGASA data, according to calculations of ™2» . The neutrino 
flux is dominated by the halo component with small admixture of extragalactic flux. The flux of 
extragalactic protons is shown by the lower curve (extr. p). 



5. SuperGZK neutrinos from Topological Defects (TDs) 

As has been first noticed by D. A. Kirzhnitz each spontaneous symmetry 
breaking in the early universe is accompanied by the phase transition. Like the phase 
transitions in liquids and solids, the cosmological phase transitions can give rise to 
topological defects (TDs), which can be in the form of surfaces (cosmic textures), lines 
(cosmic strings) and points (monopoles). In many cases TDs become unstable and 
decompose to constituent fields, superheavy gauge and Higgs bosons (X-particles), 
which then decay producing UHECR. It could happen, for example, when two seg- 
ments of ordinary string, or monopole and antimonopole touch each other, when 
electrical current in superconducting string reaches the critical value and in some 



other cases. The decays of these particles, if they heavy enough, produce particles of 
ultrahigh energies including neutrinos. 

The following TDs are of interest for UHECR and neutrinos: 
monopoles (G — > H X £7(1) symmetry breaking), ordinary strings [U(l) symmetry 
breaking) with important subclass of superconducting strings, monopoles connected 
by strings (G — > H x U{1) symmetry breaking with subsequent U(l) — > Z N symmetry 
breaking, where Zn is discrete symmetry). The important subclass of the monopole- 
string network is given by necklaces, when = Z 2 , i.e. each monopole is attached 
to two strings. We shall shortly describe the production of UHE particles by these 
TDs. 

(i) Superconducting strings. 

As was first noted by Witten"*, in a wide class of elementary particle models, strings 
behave like superconducting wires. Moving through cosmic magnetic fields, such 
strings develop electric currents. Superconducting strings produce X particles when 
the electric current in the strings reaches the critical value. Superconducting strings 
produce too small flux of UHE particles to be the sources of observed UHECR. 

(ii) Ordinary strings. 

There are several mechanisms by which ordinary strings can produce UHE particles. 

For a special choice of initial conditions, an ordinary string loop can collapse to 
a double line, releasing its total energy in the form of X-particles. However, the 
probability of this mode of collapse is extremely small, and its contribution to the 
overall flux of UHE particles is negligible. 

String loops can also produce X-particles when they self- intersect. Each intersec- 
tion, however, gives only a few particles, and the corresponding flux is very small. 

Superheavy particles with large Lorentz factors can be produced in the annihila- 
tion of cusps, when the two cusp segments overlap. The energy released in a single 
cusp event can be quite large, but again, the resulting flux of UHE particles is too 
small to account for the observations. 

It has been argued that long strings lose most of their energy not by production 
of closed loops, as it is generally believed, but by direct emission of heavy X-particles. 
If correct, this claim will change dramatically the standard picture of string evolu- 
tion. It has been also suggested that the decay products of particles produced in 
this way can explain the observed flux of UHECR However, as it is argued in 
Ref . SHI ; numerical simulations described in allow an alternative interpretation 
not connected with UHE particle production. 

(in) Network of monopoles connected by strings. 
The sequence of phase transitions 

G^H x U(l) -+H xZ N (9) 

results in the formation of monopole-string networks in which each monopole is at- 
tached to N strings. Most of the monopoles and most of the strings belong to one 



infinite network. The evolution of networks is expected to be scale-invariant with 
a characteristic distance between monopoles d = nt, where t is th e ag e of Universe 
and k = const. The production of UHE particles are considered in^*. Each string 
attached to a monopole pulls it with a force equal to the string tension, \i ~ rj 2 , 
where r] s is the symmetry breaking vev of strings. Then monopoles have a typical 
acceleration a ~ fi/m, energy E ~ \id and Lorentz factor T m ~ fid/m, where m is the 
mass of the monopole. Monopole moving with acceleration can, in principle, radiate 
gauge quanta, such as photons, gluons and weak gauge bosons, if the mass of gauge 
quantum (or the virtuality Q 2 in the case of gluon) is smaller than the monopole 
acceleration. The typical energy of radiated quanta in this case is e ~ T m a. This 
energy can be much higher than what is observed in UHECR. However, the produced 
flux (seeEHl) i s much smaller than the observed one. 
(vi) Necklaces. 

Necklaces are hybrid TDs corresponding to the case N = 2 , i.e. to the case when 
each monopole is attached to two strings. This system resembles "ordinary" cosmic 
strings, except the strings look like necklaces with monopoles playing the role of 
beads. The evolution of necklaces depends strongly on the parameter 



where m is a mass of a monopole, \x is mass per unit length of a string (tension of a 
string) and d is the average separation between monopoles and antimonopoles along 
the strings. As it is argued in Ref . , necklaces might evolve to configurations with 
r 3> 1. Monopoles and antimonopoles trapped in the necklaces inevitably annihilate 
in the end, producing first the heavy Higgs and gauge bosons (X-particles) and then 
hadrons. The rate of X-particle production can be estimated as^* 



This rate determines the rates of pion and neutrino production with energy spectrum 
calculated in Ref.^SJ. 

Restriction due to e-m cascade radiation demands the cascade energy density 
Ucas < 2 ■ 10~ 6 eV/cm 3 . The cascade energy density produced by necklaces can be 
calculated as 



where f n ~ 0.5 is a fraction of total energy release transferred to the cascade. There- 
fore, r 2 fi and the rate of X-particle production (fTTj) is limited by cascade radiation. 

The fluxes of UHE protons, photons and neutrinos from are shown in Fig. El ac- 
cording to calculations of The mass of X-particle is taken mx = 1 x 10 14 GeV. 
Neutrino flux is noticeably higher than in the case of conservative scenarios for cos- 
mogenic neutrinos and neutrinos from SHDM. 



r = m/fid, 
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Figure 5: Diffuse spectra of neutrinos, protons and photons from necklaces. The upper curve 
shows neutrino flux, the middle - proton flux and two lower curves - photon fluxes for two cases of 
absorption. The thick curve gives the sum of the proton and the higher photon flux. 



6. Mirror neutrinos 

Mirror matter can be most powerful source of superGZK neutrinos not limited by 
the usual cascade limit . 

Existence of mirror matter is based on the deep theoretical concept, which was 
introduced by Lee and Yang^^, Landau and most notably by Kobzarev, Okun 
and Pomeranchuk S3 . Particle space is a representation of the Poincare group. Since 
the space reflection x — > —x and time shift t — > t + At commute as the coordinate 
transformations, the corresponding inversion operator I s and the Hamiltonian H must 
commute, too: [I s , H] = 0. Because the parity operator P does not commute with 
H (i.e. parity is not conserved) Lee and Yang suggested that I s — P ■ R, where the 
operator R generates the mirror particle space, and thus I s transfers the left states 
of ordinary particles into right states of the mirror particles and vise versa. In fact, 
the assumption of Landau is similar: one may say that he assumed R = C. 

The mirror particles have interactions identical to the ordinary particles, but these 
two sectors interact with each other only gravitationallyl^. Gravitational interaction 
mixes the visible and mirror neutrino states, and thus causes the oscillation between 
them. 

A cosmological scenario must provide the suppression of the mirror matter and in 
particular the density of mirror photons and neutrinos at the epoch of nucleosynthe- 




Figure 6: Diffuse flux of the visible neutrinos from mirror necklaces with M x = 1 x 10 14 GeV. The 
flux is limited by observations of RICE, GLUE and FORTE. Note, that neutrino flux exceeds the 
general cascade upper limit QJ. TD curve gives the flux from the ordinary necklaces. 

sis. It can be obtained in the two-inflaton model The rolling of two inlatons to 
minimum of the potential is not synchronized, and when the mirror inflaton reaches 
minimum, the ordinary inflaton continues its rolling, inflating thus the mirror mat- 
ter produced by the mirror inflaton. While mirror matter density is suppressed, the 
mirror topological defects can strongly dominate Mirror TDs copiously pro- 

duce mirror neutrinos with extremely high energies typical for TDs, and they are 
not accompanied by any visible particles. Therefore, the upper limits on HE mirror 
neutrinos in our world do not exist. All HE mirror particles produced by mirror TDs 
are sterile for us, interacting with ordinary matter only gravitationally, and only mir- 
ror neutrinos can be efficiently converted into ordinary ones due to oscillations. The 
only (weak) upper limit comes from the resonant interaction of converted neutrinos 
with DM neutrinos: v + —> Z° . We shall obtain here this upper limit in the 
simplified case of degenerate neutrinos with common mass m v . 
The cascade energy density can be calculated as 

u cas = 27T-/ } -cr t n„ i toE$I v (E ), (13) 
hot 



where 



E 




1.81 ■ 10 



,13 




— 



2m 



V 



is the resonant neutrino energy, n Vi is the density of DM neutrinos, f to t and fhad are 
total and hadron widths of Z° decay, respectively, and 



is the effective w-cross-section in the resonance. 

Eq. (fT3j) gives the upper bound on I U (E ) which is very weak, due to factor a t n Ui t , 
as compared with that for visible neutrinos. 

The strongest limit on the fluxes of superGZK neutrinos are given nowadays by 
radio observations - . 

The mirror neutrino flux can be calculated for the case of mirror necklaces iden- 
tically to the calculations in Section 5 for ordinary necklaces, but with parameter 
r 2 /mu not being limited any more by the cascade upper bound. The probability of 
oscillation is given by P osc = 1/2, since oscillation lengths are very small in compari- 
son with the distances to TDs. The calculated neutrino fluxes for Mx = 1 x 10 14 GeV 
are shown in Fig. |H1 together with radio upper limits. The calculated flux exceeds the 
cascade upper limit for ordinary neutrino sources shown in Fig. |U1 

7. Conclusions 

SuperGZK neutrinos with energies higher than 1 x 10 20 eV can be efficiently 
searched for by future space detectors EUSO and OWL, and by radio methods. The 
neutrino-induced inclined EAS can be detected by Auger. The future detectors can 
control very large area (up to ~ 10 5 km 2 in case of EUSO) and thus they are sensitive 
to very low superGZK neutrino fluxes. The energy threshold of these methods is 
typically high, and it makes the superGZK neutrinos the main goal of the search. 

The most conservative mechanism of superGZK neutrino production is given by 
interaction of UHECR with CMB photons. One might think that the basic elements 
for UHE neutrino generation are reliably known: the beam of observed UHECR and 
the target, build by CMB photons. However, the observed flux of UHECR does 
not guarantee the detectable flux of superGZK neutrinos. A very reasonable model, 
which describes perfectly well the observed UHECR spectrum, predicts the neutrino 
flux an order of magnitude lower that of the observed UHECR flux (see upper panel 
of Fig. I2J). The detectable fluxes of superGZK neutrinos require three conditions: (i) 
the maximum acceleration energy E max > lx 10 20 eV, (ii) the cosmological evolution 
of the UHECR sources (most probably AGN) and (Hi) flat generation spectrum (e.g. 
oc E~ 2 favors the large neutrino flux). The necessary conditions (i) and (ii) imply the 
unknown astrophysics. It is especially true for (i): there are no reliable mechanisms of 
acceleration with E max ~ 10 22 — 10 23 eV, though many ideas have been put forward. 



a t = 4Sirf„G F = 1.29 • 1(T 32 cm 2 , 



(14) 



The lower panel of Fig. El presents the superGZK neutrino fluxes for the extreme 
hypothetical assumptions: very large E maK and strong evolution of the sources up to 

^max 6. 

The top-down scenarios predict naturally very high neutrino energies up to ~ 
0.1m G uT, and in some cases (monopole-string network) up to m P1 . The fluxes of 
neutrinos are also naturally high. The neutrino fluxes are rigorously constrained by 
the cascade upper limit (JJJ. The mirror neutrinos do not respect this limit, and their 
fluxes can be even larger (see Fig. [HJ) . 

The search for superGZK neutrinos is in any case is the search for a new physics, 
either for astrophysics (the new acceleration mechanisms and cosmological evolution 
of the sources, most probably AGN) or for topological defects, mirror topological 
defects and superheavy dark matter. 
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